Chaotic universe in the z = 2 Hofava-Lifshitz gravity 
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The deformed z — 2 Hofava-Lifshitz gravity with coupling constant u leads to a nonrelativistic 
"mixmaster" cosmological model. The potential of theory is given by the sum of IR and UV 
potentials in the ADM Hamiltonian formalism. It turns out that adding the UV-potential cannot 
suppress chaotic behaviors existing in the IR-potential. 
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I. INTRODUCTION 

Recently, quantum gravity at a Lifshitz point, which is 
power-counting renormalizable and hence potentially UV 
complete, was proposed by Hofava [H-Q. This theory, un- 
like string theory, is not intended to be a unified theory 
but quantum gravity. Specific cosmological implications 
of z = 3 Hofava-Lifshitz gravity with the Friedmann- 
Robertson- Walker metric based on isotropy and homo- 
geneity have recently been shown in [4-6], including ho- 
mogeneous vacuum solution with chiral primordial gravi- 
tational waves Q and nonsingular cosmological evolution 
with the big bang of standard and inflationary universe 
replaced by a matter bounce [8T4l0j. As far as the cos- 
mological solutions are concerned, there is no difference 
between z = 2 | and z = 3 | Hofava-Lifshitz gravi- 
ties because the Cotton tensor vanishes when using the 
Fricdmann-Robertson- Walker metric. Furthermore, one 
has introduced the deformed z = 3 Hofava-Lifshitz grav- 
ity to find asymptotically flat background [III (3 ■ 

On the other hand, the equations of general relativ- 
ity lead to singularities when we look at the equations 
backwards the origin of time. Especially, we concentrate 
on a temporal singularity of the solutions to the Einstein 
equations for the mixmaster model (Bianchi IX Universe) 
describing an anisotropic and homogeneous cosmology. It 
was well known that the approach to singularity shows a 
chaotic behavior. The mixmaster universe |L3l - l2C| could 
be described by a Hamiltonian dynamical system in a 
6D phase space. Belinsky, Khalatnikov, and Lifshitz 
(BKL) had conjectured that this 6D phase system could 
be well approximated by a ID discrete Gauss map that 
is known to be chaotic as one approaches the singular- 
ity [2l[. Chernoff and Barrow suggested that the mix- 
master 6D phase space could be split into the product of 
a 4D phase space and a 2D phase space having regular 
variables (H)]. Following Cornish and Levin (l7| . Lehner 
and Di Menza found that the chaos in the mixmaster 



universe is obtained for the Hamiltonian system with po- 
tential having fixed walls, which describes the curvature 
anisotropy [19]. However, it turned out that the mix- 
master chaos could be suppressed by (loop) quantum ef- 
fects mm. 

Hence it is very interesting to investigate cosmologi- 
cal application of Hofava-Lifshitz gravity in conjunction 
with the mixmaster universe based on the anisotropy 
and homogeneity because the Hofava-Lifshitz gravity is a 
strong candidate for quantum gravity. However, we will 
not make any quantum operation on the Hofava-Lifshitz 
gravity. In this work, we will make a progress on this 
direction. 

Getting an associated Hamiltonian within the ADM 
formalism 1241 of deformed z = 2 Hofava-Lifshitz grav- 
ity [H, HH EH ! we find two potentials in 6D phase space: 
IR-potential Vm from 3D curvature R and UV-potential 
Vuv from curvature square terms of R 2 and RijR lJ with 
UV coupling parameter lu. In 4D phase space, we find 
that the UV-potential cannot suppress chaotic behaviors 
existing in the IR-potential. After an extended analysis 
with movable wall, the chaotic behaviors persist in the 
6D phase space. 



II. DEFORMED z = 2 HORAVA-LIFSHITZ 
GRAVITY 

The action of the deformed z — 2 Hofava-Lifshitz grav- 
ity [H, [ll| takes the form in the (1+3)D spacetimes 



S A 



dtcPxyfgN 



.V(l-4A )i?2 _^V 
32(1 - 3A) 8 3 



(1) 



with three parameters n, n, and A. In the case of A = 1, 
the above action leads to 
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where the UV coupling parameter u> = 16/x/k is intro- 
duced to control curvature square terms [25|]. We note 
that uj is positive and thus, a negative ui is not allowed 
for the z — 2 Hofava-Lifshitz gravity. In the limit of 
uj — » oo(k 2 — > 0), 5a=i reduces to general relativity (GR) 
with the speed of light c 2 = n 2 fj, 3 /2 and Newton's con- 
stant G = k 2 /(32itc). Also, we would like to mention 
that the last line of © seems to be similar to the action 
of 3D new massive gravity when replacing /x 3 and u by 



and 



respectively [26j]. However, although a 



WwcG 

similarity between them exists, the difference is that in 
the z — 2 Hofava-Lifshitz gravity, the curvature R is a 
nonrelativistic component representing the 3D space in 
(1+3)D spacetimes, while in the 3D new massive gravity, 
the 3D curvature 3 R is a relativistic one, representing 
the (1+2)D spacetimes. Accordingly, the 3D new mas- 
sive gravity provides higher order temporal derivatives 
than second order derivatives. 

Let us introduce the metric for the mixmaster universe 
to distinguish between expansion (volume change: a) and 
anisotropy (shape change: fiij) 



ds 2 = -dt 2 



(3) 



where a 1 are the 1 forms given by 

a 1 = cos tpdO + sin tp sin 
a 2 = sin ipd9 — cos ip sin 
(T 3 = dip + cos 9d(p 



(4) 



on the three-sphere parameterized by Euler angles 
(-0, 9, cj>) with < ip < 4tt, < 9 < tt, and < 4> < 2tt. 
The shape change /3y is a 3 x 3 traceless symmetric tensor 
with det[e 2/3ij ] = 1 expressed in terms of two independent 
shape parameters f3± as 

/?n =/3+ + v / 3/3-, (3 2 2 =P+-V3p-, p 3 3 = -2/3+. (5) 

The evolution of the universe is described by giving (3± 
as function of a. Note that the closed FRW universe is 
the special case of (3± =0. 

Now we concentrate on the behavior near singularity. 
Then, the empty-space is sufficient to display the generic 
local evolution close to a singularity because the terms 
due to a matter or radiation are negligible near singular- 
ity. Using Eq. ©, the 3D curvature takes the form 



R 



12e- 2a V IR (P + ,^), 



(G) 



where the IR-potential of curvature anisotropy is given 
by 



1 r 
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2e 4/3 + cosh(4V3/3_) + e" 8/3 + 
2e~ 20 + cosh(2V3/3_) + e 4,3 + 



(7) 



The evolution of this universe is described by the motion 
of a point /3 — (j3 + , /3_) as a function of a using the time- 
dependent Lagrangian. The exponential wall picture of 
IR-potential implies that a particle (the universe) runs 
through almost free (Kasner) epochs where the potential 
could be neglected, and it is reflected at the walls, result- 
ing infinite number of oscillations. This shows that the 
system under the IR-potential behaves chaotically when 
the singularity is approached [27| . 

Before we proceed, we briefly sketch the IR-potential. 
We mention that near the point (/?+,/?_ ) = (0,0) corre- 
sponding to the global minimum, the IR-potential takes 
an approximate form of 



(8) 



As is shown in Fig. l-(a), there are three canyon lines 
located at /3_ =0 and /3_ = ±\/3/3+. We discuss 
the asymptotic structure of IR-potential. In the case 
of /3_ < 1, the IR-potential is either V IR w 2e 4 P+f3 2 _ 
if 13+ -» oo or Vir « ^e~ 8/3 + if P + -oo. For 
/3_ = 0, one has Vir if /3 + — > oo. The poten- 
tial is bounded from below and exhibits discrete Z3- 
symmetry by permuting the principal axes of rotation 
5 3 . Therefore, it has the shape of an equilateral trian- 
gle in the anisotropy space (/3+,/3_) and exponentially 
steep walls far away from the origin. A particle can only 
escape to infinity along the canyon lines where the po- 
tential has the shape shown in Fig. 2 (us — 100). The 
smallest deviation from the axial symmetry will turn the 
particle against the infinitely walls and thus, lead to a 
chaotic motion. Another useful representation of the IR- 
potential is shown in Fig. 3- (a) by drawing equipotential 
curves. They extend symmetrically between canyon lines 
at /3_ = and /3_ = ±\/3/3+, which correspond to a par- 
tially anisotropic universe with axial symmetry. The fully 
isotropic case is at the origin (0,0), where the potential 
takes the global minimum. 

The action (f5|) provides the time-dependent La- 
grangian 



C = /i 3 e 3a 



\2e- 2a V IR {P + ^„) 



16w 



V uv (p + ,I3_) 



(9) 



where the dot denotes -4r. One needs to introduce an 
emergent speed of light c in order to see the UV behav- 
iors, while for the IR behaviors, one chooses c = 1 simply. 
Here, the UV-potential of curvature square terms takes 
a complicated form 



3 



VW(A+, A-) = 4o(e 8,3 + cosh(4\/3/3_) + e 2/3 + cosh(6\/3/3_) + e~ lt)l3 + cosh(2\/3/3_)) - 40e 2/3 + cosh(2 v / 3/3-) 
+4e- 4,3+ cosh(4%/3>_) + 2e 8/3 + - 20e" 4/3 + - 42e 80 + cosh(8%/3^_) - 21e- 16,3 + 

I 



(10) 



which is a key feature of the deformed z — 2 Hofava- 
Lifshitz gravity. Note that considering the deformed 
z = 3 Hofava-Lifshitz gravity [TTJ] , one would expect to 
have a more complicated UV potential because of the 
presence of the Cotton tensor (28[. However, we have 
shown that the Cotton tensor does not change signifi- 
cantly the situations [29j |. 

In order to appreciate implications of chaotic approach 
to the deformed z — 2 Hofava-Lifshitz gravity, we have 
to calculate the Hamiltonian density by introducing three 
canonical momenta as 



P± 



dC 

dpi 



Po 



dC 

da 



12^ 3 e 3a a 



The normalized canonical Hamiltonian in 6D phase space 
takes the form 



n 6D = \(jp 2 ++p- -pI) + v a (p + ,p-,w), 



where the potential is given by 

V a ((3+,p.,uj)^e 4a (v IR 



192w 



Vuv 



(12) 



(13) 



Here Hqd — 12/x 3 e 3a H c using the canonical Hamiltonian 
H c . In this case, let us choose the parameter 12/i 3 = 1 for 
simplicity. Then, the Hamiltonian equations of motion 
are 



At =P±, P± = 



Mi 



in 



dVr 



UV 



dpi 



-Pa, p a = -4e 4a V IB + 



192w d(3± ' 

(14) 



in 6D phase space. 



III. CHAOTIC BEHAVIORS IN 4D PHASE 
SPACE 

Chernoff and Barrow showed that the mixmaster 6D 
phase space could be split into the product of a 4D phase 
space showing chaotic behavior and a 2D phase space 
showing regular behavior [l5j . Hence, we confine the 
dynamical system to a 4D phase space. Setting a = 1, 
we consider the motion of a particle (the universe) of 
coordinates (/?+,/?_) under the full potential of 



V a=1 (fi+,p-,u) -> V(P+,P-,uj) = e 4 



Vir 



Vuv 
192e 2 w 



(15) 




FIG. 1: 3D Potential shapes V(/3+, /3-,lj) with three canyon 
lines located at /3_ = and /?_ = ±\/3/3+: (a) for uj — 100, 
it corresponds to the GR mixmaster model (Vir dominates), 
(b) For uj — 0.01, it corresponds to the potential of z — 2 
Hofava-Lifshitz gravity. 



The key to chaos in mixmaster dynamics is that the 
potential has been developed "corners", that is, places 
where very special trajectories do not encounter walls. 
This leads to sensitive dependence on initial conditions, 
because a tiny change in how closely a trajectory ap- 
proaches a corner can lead to large changes in the se- 
quence of bounces off the walls. On the other hand, the 
behavior of the potential near the origin has little to do 
with chaos because the overall decrease in scale leads to 
the walls and the trajectories are moving ever farther 
outward from the origin. 

Let us describe the z = 2 Hofava-Lifshitz potential 
V(Ah A— ) w ) intensively In order to make a connection 
to the loop quantum gravity, we first mention that near 
the origin (J3+,/3-) = (0,0) which corresponds to the 
closed FRW universe, the potential V({3+, 0-,cv) takes 
approximately the form of 



Comparing (|16p with ©, the former reduces to the latter 
up to e 4 in the limit of uj — >■ oo. It turns out that adding 
the UV-potential makes just the potential well at the 
origin deeper, compared to the IR case. At this stage, 
it is curious to ask whether the inflection point at the 
origin of (/? + ,/3_) = (0,0) exists, which might show a 
change from chaotic behavior to non-chaotic behavior. 
This point may be determined by the condition of 

V"(p+,O,u>)\ 0+=o = V"(0,p-,u>)\p_ =0 = 0, (17) 



17e< 



(Pl+P 2 ). (16) 



which leads to 



2e 4 



17e 2 
2uj 



(18) 
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However, we have the negative tu c 



17 

"I? 



-0.5752. 



(19) 



This may imply that there is no inflection point which 
makes a transition from chaotic behavior to non-chaotic 
behavior. 

The asymptotic structure of the full potential is given 
as follows. For f}_ = 0, if f3 + — > 00, then 



V(/3+,P-,u) w0. 
For /3- <C 1, if P+ — > 00, one finds 

16/51 



For /3_ < 1, if P+ -> -00, then 



e- a ^G8 + ,/3_,u)« — e" 
64w 



(20) 



(21) 



(22) 



Therefore, the asymptotic structure is determined by the 
UV-potential. As is shown in Fig. l-(b), there are three 
canyon lines located at ft- = and /?_ = ± v / 3/3+- The 
potential is bounded from below and exhibits discrete 
■^-symmetry by permuting the principal axes of rotation 
S 3 . Therefore, it has the shape of an equilateral triangle 
in the anisotropy space (/?+, /3_) and exponentially steep 
walls far away from the origin. However, for total energy 
E < Vim, a particle cannot escape to infinity along the 
canyon lines where the potential has a local maximum 
V = Vim as a bump shown in Fig. 2 (ui = 0.01). For 
E < 0, the smallest deviation from the axial symmetry 
will turn the particle against the infinitely walls and thus, 
lead to a chaotic motion. Although a local maximum 
appears along canyon lines, equipotential curves in Fig. 
3-(b) are similar to Fig. 3- (a) of IR-potential for E < 
0. Hence, we expect that chaotic behavior appears in 
the z = 2 Hofava-Lifshitz gravity with small to. The 
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FIG. 2: Potential graphs V(j8+,y8_,w) with /?_ = 0: the 
long-dashed curve is for tj — 100 (GR) and the solid curve for 
cj = 0.01 with a local maximum V = Vim- 
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FIG. 3: Equipotential curves, developing triangles with cor- 
ner: (a) for lj — 100 with E — —6.5, —6.0, —5.0, —4.0, respec- 
tively, (b) for w = 0.01 with E = -15.0,-13.0,-7.0,-4.0, 
respectively. Three canyon lines are developed along corners 
at p„ = and /3_ = ±v / 3/3+. 




FIG. 4: Poincare sections for the lu = 100 case with (a) E 
-6.5, (b) E = -6.0, (c) E = -5.0, and (d) E = -4.0. 



dynamics of particle seems complicated for < E < Vi m 
and thus, we skip it. 

In general, the chaos could be defined as being such 
that (i) the periodic points of the flow associated to the 
Hamiltonian are dense, (ii) there is a transitive orbit in 
the dynamical system, and (iii) there is sensitive depen- 
dence on the initial condition. Our reduced system is 
described by the 4D Hamiltonian 

n iD ^ l -{ P l+p 2 _) + V{P + ,P-^). (23) 

It is well known that the appearance of chaotic behav- 
ior in the mixmaster dynamics is closely related to the 
appearance of "corner" potential. To make a definite con- 
nection, we choose w = 100 and u> = 0.01 for the original 
mixmaster universe of GR and the mixmaster universe of 
z = 2 Hofava-Lifshitz gravity, respectively. As is shown 
in Fig. 3, when making equipotential curves, we find that 
the corner appears for E = —5.0 and —4.0 for lu = 100 
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FIG. 5: Poincare sections for uj — 0.01 with (a) E — —15.0, 
(b) E = -13.0, (c) E = -7.0, and (d) E = -4.0. 



and for E — — 7.0 and —4.0 for ui — 0.01, which implies 
the appearance of chaos. 

Let us perform simulations of the 4D dynamics and 
represent Poincare sections, which describe the trajecto- 
ries in phase space /3+) by varying the total energy E 
or T-LiD of the system. Actually, we have performed the 
analysis for the ui = 100, 1, 0.03, and 0.01 cases. We have 
found that the chaotic behavior persists for all u> > 0. 
Figs. 4 and 5 present two typical cases, showing that 
the intersections of several computed trajectories are dis- 
placed in (p + ,f3 + ) with the plane /3_ = for different 
values of energies. In each plot, we choose initial points 
which correspond to a prescribed kinetic energy. Also we 
confirm that for w = 0.01, complicated chaotic behaviors 
appear for < E < Vi m . 

The results of Poincare sections show that consider- 
ing lower energies E — —6.5 and —6.0 for ui = 100 and 
E = —15.0 and —13.0 for ui = 0.01 within the potential 
well, the integrable behavior dominates and the inter- 
sections of trajectories represent closed curves. Impor- 
tantly, concerning higher energies E = —5.0 and —4.0 
for uj = 100 and E = -7.0 and -4.0 for u = 0.01 within 
the potential well, the closed curves are broken up gradu- 
ally and the bounded phase space fills with a chaotic sea. 
The same kinds of plots have been obtained for the other 
phase space (p_,/3_) with the plane /?+ = 0. At this 
stage, we mention that for ui — 100, the case of E = 
leads to the "corners" , where special trajectories do not 
encounter walls. However, for u> — 0.01, the case of E = 
does not lead to the "corner" because of the presence of 
a local maximum. For this, see Fig. 2. 

Finally, it confirms that the appearance of chaotic be- 
havior in the mixmaster dynamics of z — 2 Hofava- 
Lifshitz gravity is closely related to the appearance of 
corner. 



IV. CHAOTIC BEHAVIOR IN 6D PHASE 
SPACE 

Up to now, we have investigated the dynamics at fixed 
q = 1 for simplicity. This means that we have never 
explored the outer region of the potential which also de- 
termines whether the z = 2 Hofava-Lifshitz gravity is 
chaotic. We remind the reader that the true phase space 
is 6D for the vacuum universe, and thus, we have to con- 
sider a movable billiard with the potential V a (f3 + , (3-,ui) 
in Eq. (|13p because the walls are moving with time since 
the logarithm of the volume change a — i In V and its 
derivative are entering in the system. In this case, a 
and p a are regular variables as functions of time. There- 
fore, one should perform a full simulation which includes 
the dynamics of a and p a . However, this seems to be a 
formidable task and thus, we could not make a progress 
on this direction. 

In this section, instead, we will investigate a possibil- 
ity of finding chaotic behaviors by considering the small 
volume limit of a — > — oo only. To this end, it would be 
better to introduce a new time r defined by 



, dt 
= 1 V> 



„3q 



(24) 



which makes decoupling of the volume a from the shape 
At explicitly. Starting from the action @ and integrat- 
ing out the space variables, we have 



S x =i = (J- / dr 



e 3a N 
V 



6(-a 2 +/? + 2 +/3_ 2 ) 



-4a 



V 2 l2e- 2a V IR ((3+,^) 



VuviP+,/3- 



(25) 



where the prime (') denotes the derivatives with respect 
to r. Plugging N = 1 into (|25|) . we have the Lagrangian 
as 



A=i = ^ 6(-a 2 + /?; 2 + 0*) 



/2 \ 



12e 4Q ( V IR (J3+,P-) - ^ V uv^+^-) 



(26) 



The canonical momenta are given by 



/. i = %± = 12 M 3 /4, p a = ^ = -12 M V. (27) 



Then, the canonical Hamiltonian in 6D phase space is 
obtained to be 



n 



CD 



p a a' + p+f3' + + p_/3'_ - £ A =i 



-2a 



192a; 



Vuv), 



(28) 
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where we have chosen the parameter 12/i 3 = 1 for sim- 
plicity. Then, the Hamiltonian equations of motion are 
obtained as 



ol - -I 

P± = P±, P± = -e 



,,2n 



dV, 



uv 



dP± 192w df3± 

,4ai 



a' = -p a , P' a = -^e 4a V IR + —V uv . (29) 

vbuj 

We note that comparing Eqs. (|2T)1) with Eqs. (TT4"1) . there 
is no change in the Hamiltonian and its equations of mo- 
tion except replacing t by r. The evolution of a is in 
particular determined by 



, ; ,2n 



a" = Ae ia V IR - ^V uv 
9ow 



(30) 



Then, we obtain a 6D phase space consisting in the prod- 
uct of a 4D chaotic one times a 2D regular phase space 
for the a and p a variables. As the volume goes to zero 
near sing ularity (e 4Q 0, p a 0), one finds the limit 

n 6D ^^(p 2 + +p 2 _^ +K ^Had- (31) 

Hence, we note that the 6D system is not asymptotic in 
t to the previous 4D system. 

Now, we are in a position to show that the presence 
of the UV-potential does not suppress chaotic behaviors 
existing in the IR-potential. For this purpose, we have 
to introduce two velocities: particle velocity v p and wall 
velocity v w defined by 



P++P-, v w 



dr ' 



(32) 



where the wall location ft? is determined by the fact that 
the asymptotic potential K is significantly felt by the 
particle as 



P a 



2K 



e 4a-8/3 + 7e 2a-16/3+ 



12 



32uj 



(33) 



in the limit of /3+ — > — oo. On the other hand, the particle 
velocity is given by 



2H 6 D+Pl-2K. 



(34) 



In the IR-limit (uj — » oo) of Einstein gravity, the wall 
location is determined by 



"gin 



Then, the wall velocity is given by 



v IR = — 



which leads to 



dr 



,IR\ 



24p a 



(35) 



(36) 



(37) 



As a result, we find that the particle velocity is always 
greater than the wall velocity as 



,m 



2H 6D +pl-2e^V I R^\p a \>v 



IR 



(38) 



Thus, there will be an infinite number of collisions of 
the particle against the wall since it will always catch a 
wall [M HI. 

Next, let us investigate what happens in the UV-limit 
of uj — > 0. In this case, the Vjjv term dominates. The 
wall velocity takes the form 



\v» v \ = 



\Pa\ 



and the particle velocity leads to 



.uv 



\Pa\ > \V. 



UV\ 



(39) 



(40) 



in the limit of a — > —oo. This case is similar to the 
Einstein gravity. 

Finally, we could not observe a slowing down of the 
particle velocity due to the UV effects. This means that 
the chaos persists in the moving wall. 



V. DISCUSSIONS 

First of all, we point out that for uj > 0, there al- 
ways exists chaotic behavior. This contrasts to the case 
of the loop mixmaster dynamics based on loop quantum 
cosmology [23j , where the mixmaster chaos could be sup- 
pressed by loop quantum effects 22]. In the loop quan- 
tum cosmology, the effective potential at decreasing vol- 
ume labeled by "discreteness j" are significantly changed 
in the vicinity of (0,0)-isotropy point in the anisotropy 
plane — j8+. The potential at larger volumes exhibits a 
potential wall of finite height and finite extension like 
Fig. 2. As the volume is decreased, the wall moves in- 
ward and its height decreases. Progressively, the wall 
disappears completely making the potential negative ev- 
erywhere at a dimensionless volume of (2.172j) 3 / 2 in the 
Planck units. Eventually, the potential approach zero 
from below. This shows that classical reflections will stop 
after a finite amount of time, implying that classical ar- 
guments about chaos inapplicable. Once quantum effects 
are taken into account, the reflections stop just when the 
volume of a given patch is about the size of Planck vol- 
ume. 

To that end, the role of UV coupling parameter uj is 
different from the area quantum number j of the loop 
quantum gravity. In our case, time variable (related to 
the volume of V = e 3a ) as well as two physical degrees 
of anisotropy /3± are treated in the classical way without 
quantization. However, in the loop quantum framework, 
all three scale factors were quantized using the loop tech- 
niques. Hence two are quite different: the potential wells 
at the origin did not disappear for any w > in the 
z = 2 Hofava-Lifshitz gravity, while in the loop quantum 
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gravity the height of potential wall rapidly decreases until 
they disappears completely as the Planck scale is reached. 

At this stage, we compare our results with the mix- 
master universe in the generalized uncertainty princi- 
ple (GUP) [3l|. Considering a close connection between 
z = 2 Hofava-Lifshitz gravity and GUP [25|], there may 
exist a cosmological relation between them. Fortunately, 
the chaotic behavior of the Bianchi IX model was not 
tamed by GUP effects, which means that the deformed 
mixmaster universe is still a chaotic system. This is 
mainly because two physical degrees of anisotropy /3± are 
considered as deformed while the time variable is treated 
in the classical way. This supports that our approach 
without quantization is correct. 

Furthermore, it was shown that adding ( 4 i?) 2 (and 
possibly other) curvature terms to the general relati vity 
leads to the fact that the chaotic behavior is absent [32| ■ 
Hence it is very curious to see why ( 4 ii) 2 does suppress 
chaotic behavior, but —^(RijR 1 ^ — f-R 2 ) does not sup- 
press chaotic behavior. 



In conclusion, the mixmaster universe has provided an- 
other example that z = 2 Hofava-Lifshitz gravity has 
shown chaotic behavior, as other chaotic dynamics of 
string or M-theory cosmology models [30j . This may be 
because we did not quantize the Hofava-Lifshitz gravity 
and we did study its classical aspects. 
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